Exact solution of the generalized Kemmer oscillator
Zhao Zi-Long, Long Chao-Yun, Long Zheng-Wen, Xu Ting
Laboratory for Photoelectric Technology and Application, Guizhou University, Guiyang 550025, China

 

† Corresponding author. E-mail: longchaoyun66@163.com

Project supported by the National Natural Science Foundation of China (Grant No. 11565009).

Abstract

In the present work, the generalized Kemmer oscillator was introduced in one dimension. It is shown that the exact solutions of generalized Kemmer oscillator with some appropriate choices of the interactions f(x) have been obtained by using the Nikiforov–Uvarov (NU) method. Moreover, several interesting cases are discussed.

1. Introduction

In recent years, the Dirac oscillator that has been attracted a lot of attention is introduced by replacing momentum in the free particle’s Dirac equation by combination , where is the position operator, m is the particle’s mass, and is the oscillator frequency.[1] It is given the name Dirac oscillator because it reduces to the standard harmonic oscillator with a strong spin–orbit coupling in the non-relativistic limit. This system has many applications in different fields and has been studied extensively in high energy physics,[28] condensed matter physics,[911] quantum optics,[1217] and mathematical physics,[1825] etc. Especially in recent, Dutta and collaborators suggested to substitute in the one-dimensional free Dirac equation the momentum operator with and the corresponding model is called as generalized Dirac oscillator. The exact solutions of a couple of generalized Dirac oscillators with complex interactions have been obtained.[26] On the other hand, the Kemmer equation[27] is a Dirac-type equation of relativistic particle of spin-one which involves matrices obeying a different scheme of commutation rules[2830] seems to be richer than the Klein–Gordon one with respect to the introduction of interactions. Moreover, the matrices in two-body Dirac equation can be reduced to a direct sum of 10 × 10, 5 × 5, and 1 × 1 matrices realizing irreducible representations of the Kemmer algebra. Thus the DKP equation is equivalent to the two-body Dirac equation[31] and has been investigated in several active research subjects such as the quark–antiquark bound state problem.[3236] Motivated by the above basic idea of introducing generalized Dirac oscillator, in this work, we will introduce generalized Kemmer oscillator by replacing momentum with in one-dimensional free Kemmer equation and the model can be interpreted as the interaction of the anomalous magnetic moment with a general electric field. Then we will derive a complete solution of the generalized Kemmer oscillator within appropriate choices of the interactions by using the Nikiforov–Uvarov method.[3747] This paper is organized as follows. In Section 2, we introduce the generalized one-dimensional Kemmer oscillator. In Section 3, we give the exact solutions of the Kemmer equation with appropriate choices of the interaction by the Nikiforov-Uvarov method. Section 4 is devoted to the conclusion. And finally, we review the Nikiforov-Uvarov method in Appendix A.

2. Kemmer equation with generalized oscillator interaction in one dimension

The Dirac-type relativistic Kemmer equation[2732] that describes the spin-one particles is written as (in the unit system where M = 2m is the total mass of two identical spin-1/2 particles and the matrices obey the following commutation relations with In Eq. (3), is a direct product operation and , can be written respectively as where are the so-called Pauli matrices.

Now let us substitute the momentum operator with in one-dimensional free Kemmer equation (1), then the Kemmer equation with generalized oscillator interaction can be written as where the operator is chosen as with . Especially, equation (5) will reduce to Kemmer equation with general oscillator interaction for .[32,48,49] The stationary state of Eq. (5) can be described by a four-component spinor with T the transpose operator.

Putting Eq. (6) into Eq. (5), we can easily obtain four equations

From Eqs. (7), (8), (9), and (10), it is not difficult to find the following relations obeying four-component wave functions Substituting Eq. (11) into Eq. (8) will lead to the following equation In the following section, the generalized parametric NU method is used to find the solution of Eq. (12) for appropriate choices of the interactions .

3. The exact solution of Kemmer equation with appropriate choices of interactions

In this section, we concentrate our efforts in the Kemmer equation with two different appropriate choices of the interactions . The corresponding exact solutions are obtained by using the generalized parametric NU method.

3.1. The exact solution of generalized Kemmer equation with to be Coulomb potential plus linear potential

Now let us take the function to be Coulomb potential plus linear potential,[50,51] The Coulomb potential plus the linear potential has received a great deal of attention in particle physics. More correctly it is used to describe systems of quark and anti-quark bound states in the context of meson spectroscopy and it is also used to represent a radial Stark effect in hydrogen in atomic and molecular physics.[5254] In this case, equation (12) is read as where , . In order to find the possible bound state solution we introduce a new variable such as , then the equation (14) becomes Obviously, equation (15) has the same form with the basic equation of the NU method.[3747] So, it is straight-forward to gain the eigenenergy and the corresponding eigenfunction by using the NU method. By comparing Eq. (15) with Eq. (A1) in Appendix A, the following parameters can be easily given Based on the NU method, the corresponding eigenenergy and eigenfunction satisfy the following equation where is Jacobi polynomial. And the other coefficients are determined by Eq. (A4) as In Eq. (18), for the parameter , the wave function will become into the form as follows:[3747] where is the Laguerre polynomial. With the help of Eqs. (16), (17), (18), (19), and some simple mathematical manipulations, the corresponding eigenenergy and eigenfunction can be expressed respectively as follows: where N is the normalization constant. For the various values of quantum numbers n and m = 1, the square of energy spectrum of the Coulomb potential plus linear potential versus parameter a has been given in Fig. 1 (in natural units).

Fig. 1. (color online) The square of energy spectrum of the Coulomb potential plus linear potential versus a for different quantum numbers n and m = 1.

Next let us make . In the case, the interaction will become the linear potential, that is to say, the generalized Kemmer oscillator will reduce to the general Kemmer oscillator. Then, the corresponding eigenvalue and eigenfunction will change into The results are in good agreement with those given by Ref. [28].

3.2. The exact solution of generalized Kemmer equation with to be exponential-type potential

As we all know, the exponential-type potentials such as the Eckart potential,[5557] the Hulth’en-type potential,[5860] the Woods–Saxon potential,[6164] and the Morse-type potential[6568] are a typical diatomic molecular potential model which has been used to describe the interaction between two atoms in a diatomic molecule and has been widely applied in physics and chemical physics.[69] In the rest part of this paper, we will try to find the closed-form analytic solution of the Kemmer equation with the interactions to be chosen as generalized exponential potential: with , , , , and α being real parameters. If these parameters are chosen suitably, the interactions will turn into the Hulth’en-type potential, the Woods–Saxon-type potential, and the Morse-type potential. Substituting Eq. (25) into Eq. (12) gives rise to the following equation where . Then, making equation (26) can be written as where the following parameters have been defined: It is obvious that equation (27) also has the same form with the basic equation of the NU method. By comparing Eq. (27) with Eq. (A1) in Appendix A the following parameters can be easily obtained

In addition, other coefficients in the NU method are determined by Eq. (A4) in Appendix A and can be written as By using the NU method again, we can obtain the eigenfunction and eigenenergy where , , , , , and N is the new normalization constant.

Now let us consider a few special interest cases.

Fig. 2. (color online) The square of energy spectrum of generalized exponential-type potential versus α for and different quantum numbers n.
4. Conclusion

In this paper, motivated by the above basic idea of introducing Dirac oscillator, the generalized Kemmer oscillator was introduced. The generalized Kemmer oscillator has been studied when the interactions have been chosen as Coulomb potential plus linear potential and the exponential-type potentials. By using the Nikiforov–Uvarov method, the corresponding exact solutions have been found. Several special interesting cases are discussed and the results obtained agree with those obtained in previous literature. It is worth noting that apart from the cases considered in this paper, the interactions can also be taken as Mie-type potential and Kratzer–Fues potential,[70] etc.

Reference
[1] Stetsko M M 2015 J. Math. Phys. 56 012101
[2] Ito D Mori K Carriere E 1967 Nuovo Cimento 51 1119
[3] Cook P A 1971 Lett. Nuovo Cimento 1 419
[4] Moreno M Zentella A 1989 J. Phys. A: Math. Gen. 22 L821
[5] Martinez-y-Romero R P Salas-Brito A L 1992 J. Math. Phys. 33 1831
[6] Munarriz J Dominguez-Adame F Lima R P A 2012 Phys. Lett. 376 3475
[7] Grineviciute J Halderson D 2012 Phys. Rev. 85 054617
[8] Faessler A Kukulin V I Shikhalev M A 2005 Ann. Phys. 320 71
[9] Villalba V M 1994 Phys. Rev. 49 586
[10] Ju G-X Ren Z 2003 Int. J. Mod. Phys. 18 5757
[11] Quesne C Tkachuk V M 2005 J. Phys. A: Math. Gen. 38 1747
[12] Aouadi O Chargui Y Fayache M S 2016 J. Math. Phys. 57 023522
[13] Rozmej P Arvieu R 1999 J. Phys. A: Math. Gen. 32 5367
[14] Bermudez A Martin Delgado M A Solano E 2007 Phys. Rev. 76 041801
[15] Lamata L Leon J Schatz T Solano E 2007 Phys. Rev. Lett. 98 253005
[16] Sadurni E Torres J M Seligman T H 2010 J. Phys. A: Math. Theor. 43 285204
[17] Bermudez A Martin-Delgado M A Solano E 2007 Phys. Rev. Lett. 99 123602
[18] Franco-Villafane J A Sadurni E Barkhofen S Kuhl U Mortessagne F Seligman T H 2013 Phys. Rev. Lett. 111 170405
[19] Hamzavi M Eshghi M Ikhdair S M 2012 J. Math. Phys. 53 082101
[20] Lu C K Herbut I F 2011 J. Phys. A: Math. Theor. 44 295003
[21] Zarrinkamar S Rajabi A A Hassanabadi H 2010 Ann. Phys. 325 2522
[22] Chargui Y Trabelsi A Chetouani L 2010 Phys. Lett. 374 2907
[23] Bermudez A Martin-Delgado M A Luis A 2008 Phys. Rev. 77 063815
[24] Kulikov D A Tutik R S Yaroshenko A P 2007 Phys. Lett. 644 311
[25] de Castro A S Alberto P Lisboa R Malheiro M 2006 Phys. Rev. 73 054309
[26] Dutta D Panella O Roy P 2013 Ann. Phys. 331 120
[27] Kemmer N 1939 Proc. R. Soc. 173 91
[28] Boumali A 2007 Phys. Scr. 76 669
[29] Kemmer N 1938 Proc. R. Soc. 166 127
[30] Duffin R J 1939 Phys. Rep. 54 1114
[31] Bednar M Ndimubandi J Nikitin A G 1997 Can. J. Phys. 75 283
[32] Boumali A Chetouani L 2005 Phys. Lett. 346 261
[33] Merad M 2007 Int. J. Theor. Phys. 46 2105
[34] Nedjadi Y Barrett R C 1994 J. Phys. A: Math. Gen. 27 4301
[35] Chetouani L Merad M Boudjedaa T Lecheheb A 2004 Int. J. Theor. Phys. 43 1147
[36] Boutabia B Boudjedaa T 2005 Phys. Lett. 338 97
[37] Ahmadov A I Aydin C Uzun O 2014 Int. J. Mod. Phys. 29 1450002
[38] Nikiforov A F Uvarov V B 1988 Special Functions of Mathematical Physics Birkhauser Basel
[39] Hamzavi M Rajabi A A Hassanabadi H 2011 Int. J. Mod. Phys. 26 1363
[40] Hamzavi M Hassanabadi H Rajabi A A 2010 Mod. Phys. Lett. 25 2447
[41] Hamzavi M Hassanabadi H Rajabi A A 2010 Int. J. Mod. Phys. 19 2189
[42] Hassanabadi H Zarrinkamar S Rajabi A A 2013 Phys. Lett. 718 1111
[43] Rajabi A A Hamzavi M 2013 J. Theor. Appl. Phys. 7 17
[44] Rahmani S Hassanabadi H Zarrinkamar S 2014 Phys. Scr. 89 065301
[45] Kurniawan A Suparmi A Cari C 2015 Chin. Phys. 24 030302
[46] Ikot A N Obong H P Hassanabadi H 2015 Chin. Phys. Lett. 32 030201
[47] Eshghi M 2012 Chin. Phys. Lett. 29 110304
[48] Sogut K Havare A 2006 Class.Quantum Grav. 23 7129
[49] Boumali A 2004 Can. J. Phys. 82 67
[50] De Castro A S 2004 Ann. Phys. 311 170
[51] Hamzavi M Rajabi A A 2013 Chin. Phys. 22 090301
[52] Plante G Antippa A F 2005 J. Math. Phys. 46 062108
[53] Quigg C Rosner J L 1979 Phys. Rep. 56 167
[54] Chaichian M Kokerler R 1980 Ann. Phys. 124 61
[55] Dong S H Qiang W C Sun G H Bezerra V B 2007 J. Phys. A: Math. Theor. 40 10535
[56] Naderia L Hassanabadi H 2016 Eur. Phys. J. Plus. 131
[57] Soylu A Bayrak O Boztosun I 2008 J. Phys. A: Math. Theor. 41 065308
[58] Bayrak O Kocak G Boztosun I 2006 J. Phys. A: Math. Gen. 39 11521
[59] Simsek Mehmet Egrifes Harun 2004 J. Phys. A: Math. Gen. 37 4379
[60] Gonul B Ozer O Cancelik Y Kocak M 2000 Phys. Lett. 275 238
[61] Tezcan C Sever R 2009 Int. J. Theor. Phys. 48 337
[62] Berkdemir C Berkdemir A Sever R 2006 J. Phys. A: Math. Gen. 39 13455
[63] Guo J Y Fang X Z Xu F X 2002 Phys. Rev. 66 062105
[64] Ikhdair S M Sever R 2010 Int. J. Mod. Phys. 25 3941
[65] Alhaidari A D 2001 Phys. Rev. Lett. 87 210405
[66] Alhaidari A D 2002 Phys. Rev. Lett. 88 189901
[67] Bayrak G Soylu A Boztosun I 2010 J. Math. Phys. 51 112301
[68] Berkdemir C Han J 2005 Chem. Phys. Lett. 409 203
[69] Wei G F Dong S H 2010 Eur. Phys. J. 43 185
[70] Eshghi M Ikhdair S M 2014 Chin. Phys. 23 120304